Ordering and multiple phase transitions in ultra-thin nickelate superlattices 
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We interpret via advanced ab initio calculations the multiple phase transitions observed recently 
in ultra-thin LaNiOs/LaAlOs superlattices. The ground state is insulating, charge-ordered, and 
antiferromagnetic due to concurrent structural distortion and weak valency disproportionation. We 
infer distinct transitions around 50 and 110 K, respectively, from antiferromagnetic order to moment 
disorder, and from struct urally-dimerized insulator to an undistorted metallic Pauli paramagnet 
(exhibiting a cuprate-like Fermi surface). The results are in satisfactory agreement with experiment. 

PACS numbers: 71.30. +h,75.25.Dk,75.70.-i 



Incipient or actual instabilities towards collective or- 
dered states are typical of many correlated materials. 
Recent experiments (TJ [2] have investigated, via ele- 
gant nanostructuring manipulations of materials prop- 
erties, the phase transitions playing out when LaNiOa 
(LNO henceforth - in the bulk the only metallic Pauli- 
paramagnetic (PM) rare-earth nickelate) is placed in an 
intentionally perturbed environment, namely epitaxially- 
strained ultra-thin superlattice (SL) of alternating lay- 
ers of LNO and of the band insulator LaAlOa (LAO). 
For sufficiently thin LNO (2-3 layers at most), multiple 
transitions from a non-magnetic normal metal to a long- 
range-ordered magnetic, insulating, charge-ordered state 
were revealed by a crossover in conductivity tempera- 
ture (T) dependence, muon spin rotation (/iSR), spec- 
tral weight transfer in optical conductivity ^Lj, and XAS 
(x-ray absorption spectroscopy) line splitting [2]. Mag- 
netometry and /iSR strongly support long-range antifer- 
romagnetic (AF) order ^. 

The precise nature of the low-T state of the LNO/LAO 
SL and the transitions it undergoes is unclear. Here 
we address the problem from first principles studying 
a strained ultra-thin LNO/LAO SL using variational 
self-interaction-corrected local density functional theory 
(VPSIC) m [4], a parameter-free method providing an 
improved description of correlated and magnetic materi- 
als compared to semi-local approaches. So far, ab initio 
calculations have neither been able to identify an AF 
ground state as observed for these SLs, nor, a fortiori, 
to provide a general picture of the rich experimental sit- 
uation. Local (LDA) or gradient-corrected (GGA) den- 
sity functionals find neither a stable AF phase nor dis- 
torted structures. LDA/GGA+U predicts ferromagnetic 
ground states not seen in experiment, apparently irre- 
spective of U [5 . The satisfactory agreement of VP- 
SIC with experiment demonstrated below suggests an 
improved account for on-site correlations, also indicated 
by points of agreement with dynamical mean field the- 
ory (DMFT) [MS] 5 especially with the interpretation pro- 
posed in [9] during the review of the present work. 

The ground state, labeled AFD henceforth, is struc- 



turally dimerized, weakly charge-ordered, insulating, and 
AF with in-plane modulation similar to bulk rare-earth 
nickelates [10 . From calculated energies, we infer 
magnetic-ordering and metal-insulator transitions at two 
distinct critical temperatures of around 50 K and 110- 
150 K respectively. The transitions are driven by coop- 
erative structural distortion and partial valency dispro- 
portionation of Ni atoms at low T, inducing magnetic 
superexchange and Mott localization. In addition, we 
calculate the T-dependent SL conductivity within Bloch- 
Boltzmann theory [11 and discuss a possible concur- 
rent metal-insulator-transition mechanism. The results 
are consistent with the data of Ref.T, We find that the 
high-T metallic Pauli-PM phase has a Fermi surface geo- 
metrically (though not orbitally) akin to optimally-doped 
cuprates, as suggested earlier [7(112]. 

Method - Total energy, force, and bands calculations 
are performed by VPSIC [3] using the plane-wave ul- 
trasoft pseudopotential method (energy cutoff 30 Ryd, 
4x4x4 k-mesh, Gaussian smearing of 20 mRyd for the 
metal phases) in an 80- atom 2\/2x2V^x2 perovskite su- 
percell (40-atom 2^^2x^/2x4 only for the calculation of 
vertical magnetic coupling). The LNO/LAO (1+1) SL is 
simulated at the in-plane lattice constant 4.02 A, corre- 
sponding to a tensile planar strain of about 4% relative 
to the LNO bulk lattice constant. The LNO layer con- 
tains four Ni atoms to simulate AF structures. We opti- 
mize the cell length, and atomic positions via quantum 
forces [3]. The dc conductivity is calculated in Bloch- 
Boltzmann theory [11 with a relaxation time approxi- 
mation using ab initio band energies on dense k grids 
(over 2000 points). 

Structure and charge ordering - In the AFD ground 
state the nominally trivalent Nini atoms of LNO are ac- 
tually inequivalent in pairs due to a strong cooperative 
"checkerboard" dimerization of the Ni-centered octahe- 
dra (FigjTj). Ni-0 bonds alternatively expand and con- 
tract from 2 A to 2.19 A or respectively 1.83 A in a pre- 
dominantly breathing mode. Calculated octahedra rota- 
tions are minor ('^1-2°). The short and long bonds match 
those, respectively, in peroxonickel complexes with nom- 



inally tetravalent Nijv and NiO with nominally divalent 
Nin- Although the charge transfer is small, we adopt 
this labeling convention for clarity. The distortion is in- 
deed accompanied by charge transfer from Nijv to Nin, 
which we quantify by VPSIC occupations ^. The total 
transfer is 0.07 |e|, in fair agreement with 0.03 |e| esti- 
mated in [1]. This charge-ordered bond-dimerized phase 
is quite compatible with the splitting in the SL XAS spec- 
tra, analogous to insulating nickelates [2 , although we 
cannot provide a quantitative estimate of that splitting 
from our calculation. The magnetic order and insulating 
character of this state tend to confirm this conclusion, as 
discussed below. 




FIG. 1. (Color on line) Magnetic pattern of the LNO layer 
of LNO/LAO (1+1). Expanded and contracted octahedra 
around Nin and Niiv, respectively, are red (filled) and blue 
(shaded). Oxygens (not shown) sit at the shared vertexes of 
the octahedra. Dashed: supercell lattice vectors. Drawing is 
approximately to scale. 

Magnetism - The structural and valency dimerization 
are associated with the magnetic pattern in FigjT] Our 
single-LNO-layer cell has no vertical modulation by con- 
struction; the planar modulation q= (0,1/2) in the recip- 
rocal basis of the supercell (FigjT]) is the same as in bulk 
monoxides [3], and it is analogous to insulating nickelates. 
The main difference is that Nijv is entirely non-magnetic 
here, while the homologous "Ni(3-(5)" are still polarized 
in e.g. NdNiOs [10]. Indeed, Nin's [circles in the Figure] 
carry a moment /iNiii=±l-44 /i^, while Nijv's [crosses in 
the Figure] have zero moment, confirming a qualitative 
picture of Nini disproportionation into unpolarized Nijv 
^2^6^ and polarized Nijv ^2^^^. Overall, our charge-spin 
order pattern matches closely the mechanism sketched 
in Fig. lb of Ref.[13 ; on-site exchange and structural en- 
ergy gains overrule the effective on-site repulsion, quite 
screened due to the e^ states delocalization (which is, in 
turn, coherent with LNO being paramagnetic, and with 
a weakened Jahn- Teller effect). Our result clearly agrees 



with the AF long-range order suggested by magnetome- 
try and /iSR. 

The AFD phase is governed by the in-plane couplings 
Jl and Js (Fig IT]), and vertical coupling J± across the 
LAO layer. An AF Jl is expected due to superexchange 
between partially- filled Cg states. Js would be AF for 
purely x'^-y'^ Cg hopping, but as Cg states are mixed Js 
may well be FM and small. Using the energies of the 
AFD, FM, and AF-G phases in the expressions 



^r 



^F = (16 Jl + 8J5)/iLp Eg - Ep 



16J5/iNi„ 



we extract Js=-4:.l meV and a small FM 7^=0.9 meV. 
The energy difference £^A~^dF=8J±MNiii ^^ ^^^ ^-^" 
A (LAO-separated AF-stacked FM LNO planes) and 
double-FM (the same stacked FM) phases yields a tiny 
J^c^-0.05 meV, as expected due to suppressed hopping 
through LAO's Al p states very far from the Fermi en- 
ergy. The large coupling anisotropy a=JL/J±— 100 and 
the fact that Js does not contribute to the magnetic en- 
ergy of the AFD phase (see FigjT]) suggest using the Neel 
temperature T^=47rJLS'^/\n{a7p) of the 3D anisotropic 
AF Heisenberg model [M] as estimate of the critical tem- 
perature: we find Tn— 50 K, in good agreement with 40 
K experimentally. (As a only affects Tn logarithmically, 
its exact value is not essential as long as it is large.) This 
interpretation of the magnetic transition as spin order- 
disorder implicitly assumes that disordered spins fluctu- 
ate rapidly enough (say, frequency >1 MHz) above Tn so 
as not to be revealed by /iSR. 

We note in passing that FM and AF-G are also insu- 
lating and have structure, charge, and magnetic pattern 
largely similar to AFD; the main difference is that in the 
FM, /iNiiv=0-0'^ Mb, which we can neglect compared to 
/iNiii for our present estimate. 



P?- A 

■ I 

■ I 
, I 

J ; 'V 

2 




DOS ?.(.■];(." V..'Ce in 



FIG. 2. (Color on line) Right panel: Total density of states 
of the AFD (orange, shaded) and Pauli-PM (solid line) states 
(zero is the Fermi energy). Left panel: Fermi surface of the 
Pauli-PM phase in the 1x1 Brillouin zone. 

Electronic structure and metal-insulator transition - 
We now address electronic properties and discuss mech- 



anisms of metal-insulator transition. FigJ2] reports the 
total density of states (DOS), of AFD and Pauli-PM, 
right panel, and the Fermi surface of the latter (see be- 
low), left panel. The concurrent structural dimerization 
and magnetic superstructure open a 1.3-eV indirect elec- 
tronic gap. The key point is that the octahedra distortion 
is essential to obtain a gap: all undistorted phases are 
metallic and show no charge transfer. Further, only the 
Pauli-PM metal is stable among these, while the AF or 
FM states dimerize spontaneously. This suggests that the 
transition be associated to the structural dimerization, 
and that the transition temperature Tmi be identified as 
that at which the structure un-dimerizes thermally, with 
attendant gap closure. Using a stripped-down version of 
Vineyard's transition-state theory, we describe the ini- 
tially full population No=N(t=0) of distorted structural 
units (contracted and expanded Ni-octahedron pairs) as 
undergoing thermal activation out of the low-T ground 
state. The population N(t)=No exp {—Rt) is abruptly de- 
pleted, i.e. the system removes the distortion and hence 
the insulating character, for a sufficient Arrhenius activa- 
tion rate i?=z^o exp(— A£^//cbT). Since the Pauli-PM is 
the only stable undimerized state, we envisage an AFD- 
PM transition, and therefore use the AFD-PM energy 
difference AE=OAO eV per Ni pair. With a plausible ef- 
fective vibrational prefactor i^o=5 THz [15 , an activation 
rate R between 10~^ and 1 Hz (i.e., lifetimes between 280 
hrs and 1 sec) corresponds to 

Tmi = -AE/[kB{\nR-\nuo)] = 110 ^ 150 K 

in good agreement with 110 K experimentally [l]. In 
closing this section we note that the Niiv-Nin charge 
transfer is associated in optical experiments to a spectral 
weight depletion below 0.4 eV, identified as a "charge 
gap" pLj; our best shot at it is the (indeed somewhat 
larger) electronic gap, originating from the combined 
structure, charge and magnetic ordering. 

We note that in the AFD state the valence top states 
mostly project on polarized Nini, while conduction states 
do so on unpolarized Niiv- This "site-discriminated" 
gap opening is quite consistent with the "site selective 
Mott transition" proposed by DMFT calculations [9 , fur- 
ther suggesting that our method can produce, in specific 
instances, predictions matching those of sophisticated 
many-body methods. 

Fermi surface - Confirming earlier theoretical sugges- 
tions |71 112^ , the Pauli-PM phase has a single-sheet hole- 
like Fermi surface (FigJ2J left panel) centered at the 
1x1 Brillouin zone corner and analogous to optimally- 
doped cuprates. The states character is, however, mixed 
Cg rather than pure x^-y'^^ as also found by recent dy- 
namical mean field calculations [16]. The nearly two- 
dimensional pockets in Fig. [2] should give rise to quantum 
oscillation as function of inverse magnetic field, with po- 
tentially observable frequencies of about 20 kTesla (com- 
pare the 30 kTesla in e.g. metallic In [17 ). 



Transport - Since we cannot calculate the T depen- 
dence of dielectric response measured in [1 , we use Bloch- 
Boltzmann theory to calculate the dc conductivity, with 
the goal of associating the transition T with the zero of 
d<j/dT as suggested in pH. We tuned the energy depen- 
dence of the relaxation-time model [TT] to reproduce the 
T-dependence (not the value) of a in the metal phase 
using the Pauli-PM bands. The model was then fed the 
AFD ground-state bands to obtain its conductivity vs T. 
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FIG. 3. Calculated conductivity vs. T for the metallic PM 
phase (dashed) and n-doped AFD phase (dash-dot); solid line: 
linear interpolation of the metal and insulator curves. The 
value of 11 is fixed at the vertical dashed line in Fig|4] 

We conjecture that the un-dimerization metal- 
insulator transition will cause the conductivity to cross 
over smoothly from the insulator to the metal. This 
can only be assessed qualitatively in the present con- 
text. First, our method cannot describe the "dirty" metal 
phase, which exhibits an unusually low experimental con- 
ductivity; to account for this, we rescale the calculated 
metal <j by a factor 1/15, the ratio of the relaxation time 
[1 for the SL metal phase to that of a normal metal (Al). 
Second, the insulating phase's experimental conductivity 
is much higher than that of our undoped insulator at the 
relevant temperatures. We assume that this is due to 
a background impurity of unidentified origin, and thus 
calculate <j in the insulating phase for the chemical po- 
tential /i set to n-type. In FigJ3]we interpolate linearly 
the two cr's just discussed (dashed and dash-dotted lines) 
vs.T, obtaining a result (solid line) qualitatively similar 
to experiment ([1 , Fig. S8B of supporting material) and 
with dcr/dT=0 around 150 K, fairly consistent with un- 
dimerization. We note in passing that log a in the insu- 
lator phase is linear in 1/T as in pL^. Hopping behavior [2 
may be due to local disorder [18], which we cannot assess. 

Conductivity calculations as a function of /i, and 
specifically in n-type conditions, suggest a further pos- 
sible, concurrent transition mechanism. cr(/i,T) in Figj4] 
exhibits three distinct regimes as a function of increasing 
ji. Initially, a is insulator-like vs./i and T, saturating at 
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FIG. 4. (Color on line) Conductivity vs. chemical potential 
and T from 10 K, thickest blue (black) line to 290 K, thickest 
red (gray) line (20-K steps). Inset: DOS near conduction 
edge. Dashed (solid) vertical line marks intermediate regime 
(metallic crossover). Chemical potential zero at midgap. 



high T for fixed \i\ then, for \i just below the conduc- 
tion edge, G is insulator-like at low T, and crosses over at 
higher T to a normal-metal-like linear decrease; finally, at 
larger /i, a decreases with T (and grows with /i) linearly 
as in a normal metal. This behavior is related to the 
slope change of the AFD near-conduction DOS (inset of 
FigQ, with the crossover to metallic conduction occur- 
ring for /i above the DOS cusp (vertical solid line). The 
insulator cr(T) in Figjsjis obtained for /i just below the 
conduction edge (dashed vertical line in FigH, and has 
dcr/dT=0 at about 250 K. This calculation thus shows 
that an insulator-metal crossover can also be legitimately 
associated to a conduction-band-edge Fermi level pinning 
(e.g. by shallow-donor defects). Clearly, this mechanism 
would be preempted by, or at most concurrent with, the 
un-dimerization transition discussed earlier, which is ro- 
bustly rooted in the structural and magnetic energetics. 

Summary - LNO/LAO ultra-short-period SLs have a 
magnetic charge-ordered insulator ground state, making 
magnetic and metal-insulator transitions at temperatures 
we estimate in 50-70 K and 110-150 K respectively to a 
metal phase with a cuprate-like Fermi surface. This in- 
terpretation is in good agreement with available exper- 
imental data. While resulting from several cooperative 
effects, this state is basically produced by an instability of 
a " checkerboard" breathing mode (induced in turn, from 
a chemical viewpoint, by a valency instability). This is 
indirectly supported by recent experimental evidence of 
the inverse effect, i.e. the destabilization and metalliza- 
tion of insulating NdNiOs via phonon injection from the 
LAO substrate [19 . 
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